Resonant scattering by realistic impurities in graphene 
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We develop a first-principles theory of resonant impurities in graphene and show that a broad range 
of typical realistic impurities leads to the characteristic sublinear dependence of the conductivity 
on the carrier concentration. By means of density functional calculations various organic groups as 
well as ad-atoms like H absorbed to graphene are shown to create midgap states within ±0.03 eV 
around the neutrality point. A low energy tight-binding (TB) description is mapped out. Boltzmann 
transport theory as well as a numerically exact Kubo formula approach yield the conductivity of 
graphene contaminated with these realistic impurities in accordance with recent experiments. 



The mechanism determining the charge carrier mobil- 
ity of present graphene samples is being controversially 
debated. The main experimental fact requiring an ex- 
planation is that, away from the neutrality point, the 
conductivity of graphene is weakly temperature depen- 
dent and approximately proportional to the carrier con- 
centration n e Hill. This definitely requires the assump- 
tion of some long-range interactions with scattering cen- 
ters. The Coulomb interaction with charge impurities 
is an "explanation by default" Q. However, it seems 
that some experimental data cannot be explained in this 
way, especially, a relatively weak sensitivity of the elec- 
tron mobility to dielectric screening Thus, alternative 
scattering mechanisms are also discussed, such as frozen 
ripples [fl and resonant scatterers @-|l. In the first case 
the long-range character of the interactions is due to the 
long-range character of elastic deformations and in the 
second one due to divergence of the scattering length. 
New experimental data j9j seem to support the latter 
possibility. 

Theoretically, both suggestions face with serious prob- 
lems. The "ripple" mechanism requires quenching of the 
thermal bending fluctuations 0, but there are still 
no realistic scenarios of such a quenching. Resonant scat- 
tering naturally appears for vacancies |8| but they do not 
exist, in noticeable concentrations, in graphene samples 
if they are not created artificially, e.g., by irradiation |ll| . 
Adsorbates on graphene can provide resonances (quasilo- 
calized states) close enough to the neutrality point [l2i - 
15 1 but not necessarily 12, [l6j]. For impurity resonances 
some 100 meV off the neutrality point the conductiv- 
ity should display a pronounced electron-hole asymmetry 
16] which is not observed in experiments. So, it is not 
clear whether resonant impurity scattering can be the 
main limiting factor in a general case. 

In this Letter, we build a first-principles theory of elec- 
tron scattering by realistic resonant impurities, such as 
various organic molecules which are always present in ex- 
foliated graphene samples 17 , 18[ . Combining the Boltz- 
mann equation approach and a numerically exact Kubo 



formula consideration with first-principles parameters, 
we show that this class of impurities can limit electron 
transport in typical exfoliated graphene samples and ex- 
plain the experimentally observed concentration depen- 
dence of the conductivity. 

Exfoliated graphene samples are contaminated with 
long polymer chains 



17, \tt 



Most important about 
these contaminants is their possibility to form a chem- 
ical bond to carbon atoms from the graphene sheet. To 
model such a situation we carry out density functional 
theory (DFT) calculations of graphene with adsorbed 
CH3, C2H5, CH2OH (as simplest examples of different 
organic groups), as well as H and OH groups. From the 
resulting supercell band structures we derive effective in- 
teraction parameters entering a TB model and find that 
the exact chemical composition is not essential: the pa- 
rameters are very similar for all adsorbates except for the 
case of hydroxyl. This facilitates us to obtain the effect 
of the contamination on the electron conductivity. 

An atomistic description of the graphene adsorbatc 
systems is achieved by DFT calculations within the gen- 
eralized gradient approximation (GGA) [l9j] on 3 x 3-9 x 9 
graphene supercells containing one impurity. Using the 
Vienna Ab Initio Simulation Package (VASP) (20l | with 
the projector augmented wave (PAW) [ill HI basis sets, 
we obtain fully relaxed adsorption geometries and calcu- 
late the supercell band structures. 

The DFT results for CH 3 , C 2 H 5 , CH 2 OH on graphene 
are shown in Fig. [l}i and compared to H and OH adsor- 
bates. All of these impurities bind covalently to graphene 
and create a midgap state as characteristic for monova- 
lent impurities [l5| . For all adsorbates except OH the 
midgap state lies within ±0.03 eV around the neutrality 
point. As the supercell band structures for the organic 
groups and for H on graphene virtually coincide within an 
energy range of more than ±leV, it becomes clear that 
the parameters of the midgap state depend very weakly 
on the adsorbed group and, thus, can be considered as 
robust for further use in the transport theory. 

For an analytical description of these systems we start 
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Figure 1: (Color online) (a) Band structures 
of 4 x 4 graphene supercells with CH3, C2H5, 
CH2OH, H and OH adsorbates and the respective 
adsorption geometries of the CH3, C2H5, CH2OH 
(c-e) groups, (b) Comparison of the supercell 
band structure of graphene with CH3 as obtained 
from DFT to the TB models with V = 2t and on- 
site energies e d = — 0.16eV and e d — — 0.65eV. 



with a TB model of graphene, 



(1) 
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where Cj denotes the Fermi operator of an electron in the 
carbon p z orbital at site i, the sum includes all pairs 
of nearest-neighbor carbon atoms, and t ss 2.6 eV is 
the nearest-neighbor hopping parameter. In this frame- 
work, we consider a "non-interacting Anderson impu- 
rity" , adding to the JT]) the localized state, Hi mp = edd^d, 
with on-site energy ed and corresponding Fermi opera- 
tor d, which is coupled to the graphene bands by V — 
Vcld + R.c. 

To describe electron transport in pristine as well as 
doped graphene correctly, the analytical model has to 
recover the realistic system within an energy window of 
some 100 meV around the neutrality point. Applying the 
same supercell boundary conditions as in the DFT sim- 
ulations to the TB impurity model, we obtain the TB 
supercell band structures as depicted in Fig. Q]b. The 
band structure of graphene with a methyl group is well 
fitted with V rj 2t = 5.2 eV and e d « -t/16 = -0.16 eV. 

For the DFT band structures of all other neutral func- 
tional groups we find a good fit of TB with \V\ > 2t 
and \e d \ < O.lt 0.26 eV. The hybridization strength V 
being a factor 2 larger than t is in accordance with the 
hybridization for hydrogen ad- atoms from Ref. [16| and 
appears very reasonable, as the imp urity forms a er-bond 
with the host atom underneath |34| . The on-site energies 
I erf I obtained here are significantly smaller than the value 
€d = 1.7 eV used for H in Ref. [16| which will make our 
results for the transport properties qualitatively differ- 
ent. We note that the model parameters extracted here 
are converged w.r.t. the supercell size. 

The scattering of electrons caused by resonant im- 
purities is described by the T-matrix (for a review, 
see Ref. Q T(E) = =— where g (E) « 
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E-t d -V*g (E)- 

mN (E), with N (E) = ^-B(D-\E\) 

and D = \J y/Zwt ~ 6eV, is the local Green function of 
pristine graphene. Correspondingly, Nq(E) is the density 
of states (DOS) per spin and per carbon atom. The T- 
matrix exhibits a resonance at E(l 
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e d = which is the energy of the midgap state. The 
impurity model parameters obtained from DFT lead to 
resonances in an energy region of ±0.03 eV around the 
Dirac point, which proves consistency of our TB model 
with DFT. 

In the Boltzmann equation approach, the T-matrix 
can be used to estimate the conductivity a: a = 
(2e 2 /h)vFkpT, where vp is the Fermi velocity and kp 
is the Fermi wave vector. For a concentration of rii im- 
purities per carbon atom, the scattering rate reads as 
|S Sill r" 1 = (2ir/h) ni \T{E F )\ 2 N {E F ) and yields 
the conductivity 



(2e 2 /h)(2Trn l \T(E F )/D\' i 



(2) 



In the limit of resonant impurities with V — >■ 00, we 
obtain T -> -l/g (E) - [^ln|f for E < D. 
Hence, the conductivity reads in this limit as 



o a (2e 2 //i)-^ln 2 
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where n e = Ep/D 2 is the number of charge carriers per 
carbon atom. Eq. ([3]) yields the same behavior as for 
vacancies Q. In the case of the resonance shifted with 
respect to the neutrality point the consideration of Ref. 
leads to the dependence 



a oc (<?o i kp In kpR) 



(4) 



where ± corresponds to electron and hole doping, respec- 
tively, and R is the effective impurity radius. 

We now investigate to which extend realistic resonant 
impurities create sublinear behavior similar to Eqs. (|3]|2]) . 
To this end, we first estimate the conductivity according 
to Eq. @ for different types of impurities (Fig. [2]) . For 
the resonant scatterers from Fig. [T] (except for OH) the 
conductivity curves are expected to lie within the region 
bounded by the curves belonging to td = — 0.26eV and 
6d = 0.26 eV. These curves are very similar to V-shape 
experimental curves 0, 0, El and can be roughly fitted 
to the limit of Eqs. (|3|) and (j4|). The effective radius R 
resulting from Eq. ((3|) is R = D/Tivp sa 0.9A and has 
been also used in the fit according to Eq. (Q} in Fig. [21 
Experimentally, sublinear behavior similar to Eqs. 
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Figure 2: (Color online) Conductivity a in the Boltzmann 
approach as function of charge carrier concentration n e (in 
units of electrons per atom) for different impurities: Impuri- 
ties with hybridization V = 2t — 5.2 eV and on-site energies 
ed = —0.26, 0.26, and 1.7eV in concentration rii — 0.1%. 
(Curves are labelled by the corresponding ea.) Fits to the 
V — > oo limit of Eq. ((3]) with m — 0.06% (dashed) as well as 
Eq. (QJ with g = 0.02 A^ 1 (dash dotted) are shown. (Here, 
n e — Ep/D 2 corresponds to the clean graphene DOS.) 



has been observed [ll[ with effective impurity radii in 
the ranee of R = 2.3 - 2.9A. However, any estimation 
of effective radii should be considered only qualitatively, 
as D and R enter the conductivity logarithmically and a 
wide range of cut-offs lead to similar conductivity curves. 

The result for impurities with V — 2t and = 1.7 eV, 
which corresponds to H ad-atoms in the model of Ref. 
16], differs qualitatively from our results and from ex- 
perimental data which emphasizes the crucial importance 
of a careful first-principles determination of the model 
parameters. In our model and for the charge carrier 
concentration being varied within \n e \ < 0.003/C-atom= 
1.1 • 10 13 cm- 2 , impurities like CH 3 , C 2 H 5 , CH 2 OH, or H 
attached to graphene lead to a Boltzmann conductivity 
with one distinct minimum close to the neutrality point. 

At low charge carrier concentrations or high impu- 
rity concentrations, the Boltzmann approach becomes 
questionable. To understand the on-set of this param- 
eter regime and the behavior of the conductivity in this 
regime, we performed numerically exact calculations of 
the conductivity in the TB model (JTJ using the Kubo 
formula. [See [25[ .] The results for two types of resonant 
scatterers, adsorbed atoms with = —t/16, V = 2t re- 
sembling CH3 groups, and for vacancies are shown in Fig. 
[31 One can see that the Boltzmann equation is applicable 
only for impurity concentrations smaller than a few per- 
cent per site (already for 5% the difference in concentra- 
tion dependence is essential). The Boltzmann approach 
does not work near the neutrality point where quantum 



corrections are dominant [6j, [26|, |27| ■ In the range of con- 
centrations, where the Boltzmann approach is applicable 
the conductivity as a function of energy fits very well the 
dependence of Eq. ©, with q = 0.02A -1 , R = 0.6A for 
n, = 0.1%, and q = 0, R = 0.5A for n x = 0.1% with 
ftp = Ef/(Tivf) as in clean graphene. 

Close to the neutrality point the conductivity deviates 
from the Boltzmann equation result of Eq. ©. Boltz- 
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Figure 3: (Color online) Conductivity a as a function of 
charge carrier concentration n e (in units of electrons per 
atom) for different resonant impurity (ed = —t/16, V = 2t) or 
vacancy concentrations (n x ) : (a) m — n x = 0.1%, (b) 0.2%, 
(c) 1%, (d) 5%. Periodic boundary conditions are used with a 
sample containing (a) 8192 x 8192 and (b-d) 4096 x 4096 car- 
bon atoms. The carrier concentrations n e are obtained from 
the integral of the corresponding DOS depicted in Fig. (T3J). 
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Figure 4: (Color online) Density of states as a function of 
energy E for different resonant impurity (e^ = —t/16, V = 2t) 
or vacancy concentrations: ni{n x ) — 0.1%, 0.2%, 1%, 5%. 



mann theory is not capable of yielding a = 4e 2 /irh for 
clean graphene at the neutrality point 0, [26[. More- 
over, resonant impurities lead to the formation of a low 
energy impurity band (see increased DOS at low ener- 
gies in Fig. U). At impurity concentrations on the order 
of a few percent (Fig. [3] c,d) this impurity band con- 
tributes to the conductivity and can lead to a maximum 
of er in the midgap region. Moreover, the impurity band 
can host two electrons per impurity. For impurity con- 
centrations below ~ 5%, this leads to a plateau shaped 
minimum of width 2n^ (or 2n x ) in the conductivity vs. 
n e curves around the neutrality point. Analyzing the 
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plateau width in experimental data (similar to the anal- 
ysis for N2O4 acceptor states in Ref. Il3|) can, thus, yield 
an independent estimate of impurity concentrations. For 
chiral disorder 28|, 29[ corresponding to the resonant im- 
purities considered, here, as well as short range disorder 
[30L [3lj (anti)localization effects can become important 
in cases like graphane, where impurity concentrations are 
varied between a few percent and 100%. In clean mi- 
cron size samples with realistic impurity concentrations 
on the order of = 0.01% — 0.1% these effects present 
merely corrections: Upon doubling the simulation cell 
length (4096 x 4096 -> 8192 x 8192) at m = 0.1% the 
changes of the conductivity at the neutrality point are 
below 10%. 

Electron scattering in bilayer graphene has been 
proven to differ essentially from the single layer case in 
Ref. l32f For a scattering potential with radius much 
smaller than the de Broglie wavelength of electrons, the 
phase shift of s-wave scattering So tends to a constant as 
k — > 0. Therefore, within the limit of applicability of the 
Boltzmann equation, the conductivity of a bilayer should 
be just linear in n e , instead of sublinear dependence (U) 
for the single layer. The difference is that in the single 
layer, due to vanishing DOS at the Dirac point, the scat- 
tering disappears at small wave vectors as So(k) oc ln ], R 
(with In 2 kR on the order of 10 for typical amounts of 
doping) for resonant and as So(k) oc kR for the nonres- 
onant impurities. Contrary, in the bilayer there are no 
restrictions on the strength of the scattering and even 
the unitary limit (8n = 7r/2) can be reached at k = 0. As 
follows from Ref. |32j, a cylindric potential well of radius 
R, leads to So = tt/2 if T^Tjm = 0> where q is the wave 
vector inside the well, Jo and Iq are the Bessel functions 
of real and imaginary arguments, respectively. Thus, an 
assumption that resonant scattering is the main limiting 
factor for electron mobility in exfoliated graphene leads 
to the prediction that the dependence of cr(n e ) should 
be essentially different for the cases of bilayer and sin- 
gle layer, that is, linear and sublinear, respectively. This 
agrees with the experimental results [331 ]. 

In summary, we have demonstrated that realistic im- 
purities in graphene frequently cause quasilocal peaks 
nearby the neutrality point. In particular, for various 
organic groups the formation of a carbon-carbon bond 
results in the appearance of midgap (resonant) states 
within ±0.03 eV around the neutrality point. They can 
be described as Anderson impurities with the hybridiza- 
tion parameter of about It and on-site energies on the 
order of |e<j| < t/10. The resonant scattering model 
with these parameters describes satisfactory experimen- 
tal data on the concentration dependence of charge car- 
rier mobility for graphene. 
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